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We present an accurate nucleon-nucleon (NN) potential based upon chi- 
ral effective Lagrangians. The model includes one- and two-pion exchange 
contributions up to chiral order three. We show that a quantitative fit of the 
NN D-w&ve phase shifts requires contact terms (which represent the short 
■ range force) of order four. Within this framework, the AW phase shifts below 

300 MeV lab. energy and the properties of the deuteron are reproduced with 
high-precision. This chiral NN potential represents a reliable starting point 
for testing the chiral effective field theory approach in exact few-nucleon and 
microscopic nuclear many-body calculations. An important implication of the 
present work is that the chiral 2n exchange at order four is of crucial interest 
for future chiral NN potential development. 
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One of the most fundamental problems of nuclear physics is to derive the force between 
two nucleons from first principles. A great obstacle for the solution of this problem has been 
the fact that the fundamental theory of strong interaction, QCD, is nonperturbative in the 
low-energy regime characteristic for nuclear physics. The way out of this dilemma is the 
effective field theory concept which recognizes different energy scales in nature. Below the 
chiral symmetry breaking scale, A x « 1 GeV, the appropriate degrees of freedom are pions 
and nucleons interacting via a force that is governed by the symmetries of QCD, particularly, 
! (broken) chiral symmetry. 

The derivation of the nuclear force from chiral effective field theory was initiated by 
Weinberg |L| and pioneered by Ordonez |2| and van Kolck []3],(4| . Subsequently, many groups 
^ ■ got involved in the subject PTTOl. As a result, efficient methods for deriving the nuclear 
force from chiral Lagrangians have emerged. Also, the quantitative nature of the chiral AW 
potential has improved fi~0 |. Nevertheless, even the currently 'best' chiral AW potentials are 



too inaccurate to serve as a reliable input for exact few-nucleon calculations or miscroscopic 
nuclear many-body theory. 

The time has come to put the chiral approach to a real test in microscopic nuclear 
structure physics. Conclusive results can, however, be produced only with a 100% quanti- 
tative AW potential based upon chiral Lagrangians. For this reason, we have embarked on 
a program to develop a AW potential that is based upon chiral effective field theory and 
reproduces the AW data with about that same quality as the high-precision AW potentials 
constructed in the 1990's |TT|-|lj|. 



*On leave from University of Salamanca, Spain. E-mail: dentem@uidaho.edu 
^E-mail: machleid@uidaho.edu 
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FIG. 1. The most important irreducible one- and two-pion exchange contributions to the NN 
interaction up to order Q 3 . Vertices denoted by small dots are from C^L, while large dots refer to 
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Starting point for the derivation of the NN interaction is an effective chiral ttN La- 

(1) 



jrangian which is given by a series of terms of increasing chiral dimension fl5 
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where the superscript refers to the number of derivatives or pion mass insertions (chiral 
dimension) and the ellipsis denotes terms of chiral order four or higher. 

We will apply the heavy baryon (HB) formulation of chiral perturbation theory |TJ| in 
which the relativistic Lagrangian is subjected to an expansion in terms of powers of 1/Mjy 
(kind of a nonrelativistic expansion), the lowest order of which is 
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where we use the notation of Ref. ||16|| . For the parameters that occur in the leading order 
Lagrangian, we apply M N = 938.919 MeV, = 138.04 MeV, f m = 92.4 MeV, and g A = 
QttNN fir/M N = 1.29, which is equivalent to gi NN /^ = 13.67. 

The HB projected Lagrangian at order two is most conveniently broken up into two 
pieces, 
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FIG. 2. D-wave phase shifts of NN scattering. The predictions by the chiral model displayed 



in Fig. 1 are shown by the solid line and the ones by the Bonn tt + 2tt model [18| by the dashed 



curve. The dotted line is OPE. Solid dots represent the Nijmegen multi-energy np analysis [19] 



and open circles the VPI/GWU analysis [20|. 
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Note that C^n 

fix is created entirely from the HB expansion of the relativistic C nN and 

thus has no free parameters ("fixed"), while ct is dominated by ttN contact terms 
proportional to the c, parameters, besides some small 1/Mn corrections. The parameters q 
are known as low-enery constants (LECs) and must be determined empirically from fits to 
ttN data. We use the values determined by Biittiker and MeiBner |17J which are (in units 
of GeV" 1 ) d = -0.81, c 3 = -4.70, and c 4 = 3.40 (c 2 will not be needed). 

The ixN Lagrangian is the crucial ingredient for the evaluation of the pion-exchange 
contributions to the NN interaction. Since we are dealing here with a low-energy effective 
theory, it is appropriate to analyze the contributions in terms of powers of small momenta: 
(Q/A x ) u , where Q is a generic momentum or a pion mass and A x ^ 1 GeV is the chiral 
symmetry breaking scale. This procedure has become known as power counting. For the 
pion-exchange diagrams relevant to our problem, the power v of a diagram is determined by 
the simple formula 
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where I denotes the number of loops in the diagram, dj the number of derivatives involved 
in vertex j, and the sum runs over all the vertices of the diagram. 

The most important irreducible one-pion exchange (OPE) and two-pion exchange (TPE) 
contributions to the NN interaction up to order Q 3 are shown in Fig. 1; they have been 
evaluated by Kaiser et al. using covariant perturbation theory and dimensional regular- 
ization. One- and two-pion exchanges are known to describe NN scattering in peripheral 
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partial waves. In G and higher partial waves (orbital angular momentum L > 4), there is 
good agreement between the chiral and conventional |18| 2-k model as well as the empirical 



phase shifts [fT9| , p0f . The agreement deteriorates when proceeding to lower L. While in F 
waves the agreement between the chiral model and the empirical phase shifts is still fair, 
substantial discrepancies emerge in D waves, Fig. 2, where the chiral 2tt exchange is far too 
attractive — a fact that has been noticed before f|[7|]. 

To control the D (and lower) partial waves, we need (repulsive) short-range contributions. 



In the conventional meson model JT8| , these are created by the exchange of heavy mesons 



(notably, the u meson). In chiral perturbation theory (xPT), heavy mesons have no place 
and the short-range force is parametrized in terms of contact potentials, which are organized 
by powers of Q. If Q is, e. g., a momentum transfer, i. e., Q = p ' — p, where p and p ' are 
the CM nucleon momenta before and after scattering, respectively, and 9 is the scattering 
angle, then, for even v, 

<T~(cos#) m with rn<^. (8) 
Partial-wave decomposition for orbital-angular momentum L yields, 

[ +1 Q u P L (cos9)dcos9 /0 for L<-, (9) 
J—i 2 

where Pl is a Legendre polynominal. The conclusion is that for non-vanishing contributions 
in D waves (L = 2), v = 4 is required. This one important message that we like to convey 
in this letter. Based upon invariance considerations, there are a total of 24 contact terms up 
to order Q 4 , which we all include in our model. The parameters of these terms have to be 
natural, but are otherwise unconstrained and, thus, represent essentially free parameters. 

To describe NN scattering, we start from the Bethe-Salpeter (BS) equation |21] which 
reads in operator notation 

t = v + vgt (io) 

with T the invariant amplitude for the two-nucleon scattering process, V the sum of all 
connected two-particle irreducible diagrams, and Q the relativistic two-nucleon propagator. 
The BS equation is equivalent to a set of two equations: 

T = V + VgT (11) 
V = V + V {G -g)V (12) 
« V + Vope (G - g)V OP E (13) 

where the last line states the approximation we are using, exhibiting the way we treat the 
27r box diagram. This treatment avoids double counting when V is iterated in the scat- 
tering equation and is also consistent with the calculations of Ref. ||. For the relativistic 
three-dimensional propagator g, we choose the one proposed by Blankenbecler and Sugar 
f22fl (BbS) which has the great practical advantage that the OPE (and the entire poten- 
tial) becomes energy- independent. Thus, we do not need the rather elaborate formalism of 
unitary transformations |10| to generate energy-independence of the potential. 
Our full chiral NN potential V is defined by 
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FIG. 3. Phase shifts for J < 2. The solid line is the result from our chiral NN potential, while 
the dotted and dashed lines are the predictions by two chiral models developed by Epelbaum et 
al. (NLO and NNLO, respectively). Solid dots and open circles represent phase shift analyses 
explained in the caption of Fig. 2. 
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where the first term on the r.h.s. is given by Eq. (|13|) with V containing essentially the 
diagrams of Fig. 1. This potential satisfies the relativistic BbS equation, Eq. (JTTJ) . If we 
define now, 
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with E p = yMj + p 2 , then V satisfies the usual, nonrelativistic Lippmann-Schwinger (LS) 
equation. 

Iteration of V in the LS equation requires cutting V off for high momenta to avoid 
infinities. Therefore, we regularize V in the following way: 



V(p',p) i — > V(p',p) e 

^V(p',p) < 1 
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where the last equation is to indicate that the exponential cutoff does not affect the order to 
which we are calculating, but introduces contributions beyond that order. For the contact 
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TABLE I. Two- and three- nucleon low-energy data. 





Idaho-A 


Idaho-U 


f ' T \ T ) [1/1 

CD-Bonn |14 


nrin [M oil 

AV18 |12|] 


Empirical^ 


Low-energy np scattering 












1 5o scattering length (fm) 


-23.75 


-23.75 


-23.74 


-23.73 


-23.74(2) 


1 Sq effective range (fm) 


2.70 


2.70 


2.67 


2.70 


2.77(5) 


3 Si scattering length (fm) 


5.417 


5.417 


5.420 


5.419 


5.419(7) 


3 Si effective range (fm) 


1.750 


1.750 


1.751 


1.753 


1.753(8) 


Deuteron properties 












Binding energy (MeV) 


2.224575 


2.224575 


2.224575 


2.224575 


2.224575(9) 


Asympt. S state (fm -1 / 2 ) 


0.8846 


0.8846 


0.8846 


0.8850 


0.8846(9) 


Asympt. D/S state 


0.0256 


0.0255 


0.0256 


0.0250 


0.0256(4) 


Deuteron radius (fm) 


1.9756 c 


1.9758 c 


1.970 c 


1.971 c 


1.9754(9) d 


Quadrupole moment (fm 2 ) 


0.281 6 


0.284 e 


0.280 e 


0.280 e 


0.2859(3) 


.D-state probability (%) 


4.17 


4.94 


4.85 


5.76 




Triton binding (MeV) 


8.14 


8.02 


8.00 


7.62 


8.48 



"Chiral NN potential of the present work. 

For references concerning the empirical data, see Tables XIV and XVIII of Ref. ||. 
c With meson-exchange current (MEC) and relativistic corrections |23|| . 
^Reference @. 

including MEC and relativistic corrections in the amount of 0.010 fm 2 f25| . 



terms, we use partial wave dependent cutoff parameters A m 0.4 — 0.5 GeV which brings 
the total number of parameters in our chiral AW model up to 46. At first glance, this may 
appear to be a large number. Note, however, that the Nijmegen phase shift analysis [H5| 
and the high-precision potentials [|ld|-|i"4|j developed in the 1990's carry between 40 and 50 
parameters. Thus, the number of parameters needed for a quantitative chiral AW model is 
just about the same as for meson models. Since the chiral model has less predictive power 
than the meson model this should not be unexpected. 

In Fig. 3, we show the phase shifts of neutron-proton (np) scattering for lab. energies 
below 300 MeV and partial waves with J < 2. The solid line represents the result from the 
chiral AW potential developed in the present work. The reproduction of the empirical phase 
shifts by our model is excellent. For comparison, we also show the phase shift predictions by 
two chiral models recently developed by Epelbaum et al. [|TU| (dotted and dashed curves in 
Fig. 3). Also the effective range parameters in S waves agree accurately with the empirical 
values, as well as the deuteron parameters (see Table I). We note that our present chiral 
potential is charge-independent and adjusted to the np data. 

Due to the very quantitative nature of this new chiral AW potential, it represents a 
reliable and promising starting point for exact few-body calculations and microscopic nuclear 
many-body theory. 

A crucial finding of our investigation is that contact terms of order four are required for 
a quantitative AW model. The basic ideas of ;\TT may then suggest that — for reasons of 
consistency — the chiral 2n exchange contribution should also be included up to order four. 



Therefore, an implication of the present work is that the chiral 2n exchange at order four |26 
will be important for further chiral AW potential development. 
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